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Abstract 



We consider a holographic dual of hydrodynamics of M = 4 SYM plasma that undergoes 
non-isotropic three-dimensional expansion relevant to RHIC fireball. Our model is a nat- 
ural extension of the Bjorken's one-dimensional expansion, and it describes an elliptic flow 
whose V2 and eccentricity are given in terms of anisotropy parameters. Holographic renor- 
malization shows that absence of conformal anomaly in the SYM theory constrains our 
local rest frame to be a Kasner spacetime. We show that the Kasner spcetime provides 
a simple description of the anisotropically expanding fluid within a well-controled ap- 
proximation. We also find that the dual geometry determines some of the hydrodynamic 
quantities in terms of the initial condition and the fundamental constants. 
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1 Introduction 



There has been much hope that one might be able to use AdS/CFT [T] to describe the real systems 
after certain amount of deformations. In fact it has been suggested that the fireball in Relativistic 
Heavy Ion Collision (RHIC) can be explained from dual gravity point of view [21 El El E], since the 
quark-gluon plasma (QGP) created there are in the strong coupling region [HE]- Although the YM 
theory described by the standard AdS / CFT is large- N c N = A SYM theory, there are many attempts 
to construct models closer to QCD [8|. SUSY is not very relevant in the finite temperature context 
because it is broken completely. 

Since the RHIC fireball is expanding, we need to understand AdS/CFT in time dependent situa- 
tions. Recently, Janik and Peschanski [91 [10] discussed this problem in non- viscous cases. They use 
the conservation law and conformal invariance together with the holographic renormalization [T2] 
to express the bulk geometry from the given boundary data. As a result, the bulk geometry repro- 
duces the basic features of Bjorken theory [13]. Theses results were generalized to the cases where 
shear viscosity is included [H] . Indeed, it had been pointed out that inclusion of shear viscosity is 
very important in the analysis of real RHIC physics since it plays an essential role in the elliptic 
flow (see for example, [151 EH1 El HH1 US GO])- The shear viscosity at the strong coupling limit was 
calculated for the M = 4 SYM systems in [21] using AdS/CFT, which fits the perfect fluidity of 
RHIC QGP (see [3 [H [22] and the references therein). 

In this paper, we consider a holographic dual of strongly interacting N = 4 large- N c SYM fluid 
with non- isotropic three-dimensional expansion which is relevant to "Little Bang" of RHIC. We will 
first make a simplest generalization of the Bjorken's one-dimensional expansion to three-dimensional 
cases. The resulting local rest frame of the fluid is described by the Kasner metric whose extreme 
limit reproduces the Bjorken hydrodynamics. Interestingly, the hydrodynamic equation and the 
equation of state will be shown to be independent of the non-isotropy parameters. 

Further more, using holographic renormalization, we will establish the holographic dual geometry 
of the anisotropically expanding fluid in the Kasner spacetime and discuss its physical consequences 
by extending the work of Q3] . Interestingly, the gravity dual carries more information than the hy- 
drodynamics which is our input. The integration constants coming from the hydrodynamic equation 
cannot be determined by the hydrodynamics itself. However, we will show that such quantities can 
be determined by considering the dual geometry. Therefore it provides a method to extract useful 
information on the macroscopic properties of the expanding fluid in terms of microscopic data. 

Kasner spacetime is a curved spacetime in general. However, we find that the Kasner spacetime 
can be a well-controled approximation of a local rest frame of an anisotropically expanding fluid on 
Minkowski spacetime. Therefore, the Kasner spacetime provides a useful framework for analysing 
the three-dimensional expansion because of the simplification of the hydrodynamic equations. 

The organization of the paper is the following. In Section [21 we introduce Kasner spacetime as 
a local rest frame of a fluid under three-dimensional anisotropic expansion. We also establish the 
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hydrodynamics based on that local rest frame. Section [3] gives analyses in the gravity dual in the 
late time regime. We extend the results of [£l Q3] to the case of our interest: the three-dimensional 
anisotropic expansion with shear viscosity. We show that the dual geometry determines some of the 
hydrodynamic quantities in terms of the initial condition and the fundamental constants. In Section 
H]we show that the hydrodynamics on the Kasner spacetime (and hence its gravity dual) describes an 
elliptic flow on the flat spacetime within a well-controlled approximation. We analyze the properties 
of the flow. We conclude in the final section. The definition of the approximation and its justification 
are given in appendix. 

2 Hydrodynamics in Kasner space and anisotropic expansion 

The remarkable success of Bjorken's hydrodynamics needs to be extended for more realistic cases of 
three-dimensional expansions. The Bjorken's basic assumption is existence of the so-called "central 
rapidity region" (CPR). The local rest frame of the fluid can be given by propertime(r)-rapidity(y), 
whose relationship with the cartesian coordinate is (X°, X , X 2 , X 3 ) = (r coshy, r sinhy, X 2 , X 3 ). 
We have chosen the collision axis to be in the x 1 direction. The Minkowski metric in this coordinate 
has the form 



Our starting point is an observation that the above metric describes a one-dimensional Hubble 
expansion, i.e., an expansion of the universe instead of that of the fluid, in which the rapidity plays 
the role of the co-moving coordinate. The hydrodynamic equations can be derived from the covariant 
conservation of the energy-momentum tensor with the above metric. 

The real expansion in RHIC is not one-dimensional but a three-dimensional one, although the 
dominant flow is along the collision axis. The idea is that since the one-dimensional approximation 
of RHIC fireball expansion is seen as a Hubble flow, a three-dimensional Hubble flow may describe 
the RHIC plasma better. Let us begin with an ansatz for a local rest frame given by 



Here x l ,i = 1,2,3 denote the co-moving coordinates. The a, b and c are arbitrary constants for a 
moment. However, as we will show in Section [3l they satisfy 



if we impose conformal invariance of the fluid. Under the above conditions, the metric is called 
Kasner metric describing a homogeneous but anisotropic expansion of the Universe. In short, we 
identify the "Little Bang" in "SYM version of RHIC" with a Big Bang with the homogeneous but 
anisotropic expansion described by the Kasner metric. We use this metric in the late time regime 
because we describe only the late time evolution of the fluid where hydrodynamics is valid. Therefore 



ds 2 = -dr 2 + r 2 dy 2 + (dX 2 ) 2 + {dX 3 ) 2 . 






a + b + c=l, a 2 + b 2 + c 2 = l, 
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the initial singularity of the metric is not a relevant feature for us. Since we take x 1 as the longitudinal 
direction, a realistic set up is to choose a ~ 1 and b, c ~ (see section 

We first establish the hydrodynamics in the Kasner metricQ We assume that the expansion is 
anisotropic but homogeneous, and the physical quantities depend only on r. Using the fact that the 
energy-momentum tensor is diagonal on the local rest frame, and using the symmetry in transverse 
coordinates, we can write 

r^ = diag(-p,/i,/ 2l / 3 ), (2.4) 
where p is the energy density of the fluid. By use of the conservation law \7 ^T^ v = we get 

p+p/T + Y, a ifi/ T = > ( 2 - 5 ) 

i 

and from the conformal invar iance T£ = 0, we get 

-P + J2fi = °- ( 2 - 6 ) 

To get the above results, we have used the following non-zero T^'s: 

Tl x =ar 2a -\ Tl 2 = br 2b -\ 1% = cr 2c -\ F\ T = a/r, rf T = 6/r, T\ T = c/r. (2.7) 

On the other hand, the energy-momentum tensor in the framework of relativistic hydrodynamics 
is known to be 

r = (p + P)uY|Pf It^, (2.8) 

where P is the pressure of the fluid, = (7, ^v) is the four-velocity field in terms of the local fluid 
velocity v, and t^ u is the dissipative term. In a frame where the energy three-flux vanishes, t^ u is 
given in terms of the bulk viscosity £ and the shear viscosity r] by 

= -ri(&» x Vxu v + A^Vau^ - \ V x u x ) - £ A^ V x u x , (2.9) 

under the assumption that r^ v is of first order in gradients. We have defined the three-frame projec- 
tion tensor as = g^ v + u^u v . For the conformal invariance, we set £ = 0. Notice that the bulk 
viscosity in the realistic RHIC setup is also negligible. (See for example, Ref. [T6].) 

The four-velocity of the fluid at any point in the local rest frame is = (1,0,0,0), and this 
makes the energy-momentum tensor diagonal. Using = diag(0, r~ 2a , t~ 26 , t~ 2c ), V\u u = = 
diag(0, a/r, b/r, c/t) and V u u u = (a + b + c)/t, we get the mixed energy-momentum tensor: 



/ -p \ 



P-f(3a-l)2 

P- f(36- 1)2 

V P-§(3c-l)2 J 



(2.10) 



1 We use the late time approximation, which is employed in |14) . where the macroscopic quantities are assumed to 
evolve sufficiently slowly. 
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By identifying (|2.10|) with (|2.4|) . we obtain 



/i =?-(<*- (2.11) 
o r 



p+ 1 p+E-^-(^-^)-)= ' ( 2 - 12 ) 



where o« = a,b,c for z = 1,2,3 respectively. Inserting these into (|2.5p and (|2.6p . we get 

~p + E — [p- q )— 

and 

- p + 3p - 27? ( Ei< ^~ 1 ) = 0. (2.13) 
Using the second equation, we may write the first equation in terms of the energy density as 

As we will discuss in Section O if we impose the conformal invariance, we get the conditions 

E«i = i, E« 2 = 1 ' ( 2 - 15 ) 

under which the equation of state and the conservation law become: 

P 

P = 3' 



3Z +oZ = o^2' ( 2 - 16 ) 



4/9 4 7/ 

cfr 3 r 3 r 2 

Remarkably, the equations for fluid dynamics are completely independent of the parameters a, b, c in 
this case. In fact the dynamical law should not depend on the initial conditions (a, b and c) that 
reflects the initial collision geometry. In this respect, (|2.15p leads us a satisfactory consequence. 

Notice that both of p and n depend on the proper time r in general. Let's assume that the shear 
viscosity evolves by rj = tjq/t 13 , where 770 is a positive constant. The solution of (|2.16p is then given 
by 

= ^ + T7^7W 1/3), ( 2 -17) 

- P ° +^ (for P= 1/3), 



T 4/3 3 r 4/3 

where po is a positive constant. For (3 < 1/3 case, the viscous corrections in the hydrodynamic 
quantities become dominant in the late time, which invalidates the hydrodynamic description. If 
(3 > 1/3, the shear viscosity term is subleading in the late time as we expect. Therefore we will 
consider only the latter case from now on. 

The proper-time dependence of the temperature T can be read off by assuming the Stefan- 
Boltzmann's law p oc T 4 : 
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In the static finite temperature system of strongly coupled N = 4 SYM theory, it is known that 
rj oc T 3 [21]. Let us assume that the result is valid in the slowly varying non-static cases. Then we 
can set (3 = 1 hence 

r, = ^. (2.19) 
r 

We know p ~ T 4 and rj ~ T 3 can be consistent only if there is an additional term in (|2. 19[) , but the 
correction term is negligible in our case. The temperature behavior is then given by 

T = T Q ( * -^V-V (2.20) 



-V3 2 pq t 

We can evaluate the entropy change in the presence of shear viscosity by using hydrodynamics. 
The conservation of energy-momentum tensor can be rewritten as 

dr + dr ^ ~ 3 r2 ' [Z - ZL) 
where yfg = r is the volume element in the co-moving coordinate. By integrating over the unit 
volume in the co-moving coordinate, and using the thermodynamic relation between the entropy and 
energy- work, 

dE + PdV = TdS, (2.22) 

(I2.21|) can be written as 



d(y/gs) _Arj^g 
1 ~ 3~ [Z - 26) 



where s denotes the entropy density and ^fgs = S is the entropy per unit co-moving volume. Notice 
that in the absence of viscosity, the entropy per unit co-moving volume is constant. Now, the entropy 
per unit co-moving volume has the time dependence: 

S(r) = S + 4 -[ T dr^ 



3 Jo t 2 T 

= Soo -2^r^ + ..., (2.24) 

where Soo = S(oo). The dissipation creates entropy but its rate slows down with time. Notice that 
all these arguments are completely in parallel with the case of the one-dimensional expansion. 

The hydrodynamics does not calculate the constants and To in terms of the initial condition. 
It is important to point out that by embedding the hydrodynamics into AdS/CFT, we can determine 
these parameters in terms of the initial condition p$ and the fundamental constants of the theory. 



3 Holographic dual of anisotropic expansion 

In this section, we will find a five-dimensional metric that is dual to the hydrodynamic description of 
the YM fluid in the previous section. Some of the parameters of hydrodynamics will be determined 
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as a consequence. The basic strategy is to use the Einstein's equation together with the boundary 
condition given by the four-dimensional energy-momentum tensor [T2"l [9] . We consider general 
asymptotically AdS metrics in the Feffer man- Graham coordinate: 

, 2 2 9^dx> 1 dx' / + dz 2 

as =r -z g ' I 3 - 1 ) 

where = (t, x 1 , x 2 , x 3 ) in our case, tq = (Ang s N c a' 2 ) l / 4: is the length scale given by the string 
coupling g s and the number of the colors N c . The four-dimensional metric g^ u is expanded with 
respect to z in the following form I12j : 

9flu (T, z) = g$(r) + z 2 g$(r) + z 4 g$(r) + z & g$(r) + ■■■. (3.2) 

Here gffi is the physical four-dimensional metric for the gauge theory on the boundary, which is 
given by (|2.2|) in the present case. The fl^'s depend only on r because our physical quantities are 
assumed to depend only on r. 

gffi is related to the conformal anomaly of the YM theory in the following way 

W = ~vtcl [V**®) 2 ~ T^) 2 ] > (3-3) 
where G5 is the 5d Newton's constant given by G5 = 87r 3 a' 4 g 2 /rQ so that 

4ttG 5 2vr 2 



" N 2 



(3.4) 



in our notation. Since we are dealing with M = 4 SYM theory on R ' 3 , we should require the 
conformal invariance: (Ttf) = 0. The most natural choice is given by 

9$ = 0. (3.5) 

This is equivalent to the Ricci flat condition for the four-dimensional metric: 

iV = 0, (3.6) 



{ $ = ~[R»v-~R9m. (3.7) 



through the relationship 



Since 

^oo = (E^-E a2 )/ r2 ' Ru = a l C£a j -l)T 2a >- 2 , (3.8) 



the Ricci flat condition (|3.6|) gives the Kasner condition 

One should notice that we are not solving the four-dimensional Einstein's equation in the presence 
of T^y. The Kasner metric is not a consequence of the gravitational effect of T^, but an effective 



6 



description of the spacetime expansion satisfying the conformal invariance of strong gauge theory 
interaction. 

We can identify the first non-trivial data in J32D, g$, with the energy-momentum tensor at the 
boundary [TT] : 

4vrG 5 



9 (4) 



(3.9) 



in our notation. For the time being, we set ^ttGz, = 1 and rg = 1. The higher-order terms in (|3.2p 
are determined by solving the Einstein's equation with the negative cosmological constant A = —6 
[D] (see also [9]): 

1 



Rmn — t;GmnR — 6Gmn = 0, 



(3.10) 



where the metric and the curvature tensor are the five-dimensional ones of (|3.ip . g^v is described 
by 9^v 2 \g^v j • • • ,9$ through the Einstein's equation in five dimensions. In other words, we 
can obtain the higher-order terms in (|3.2p recursively by starting with the initial data gffu (~Kasner) 
and g$ (~ T^ u ). 

Let us come back to our main interest to obtain the bulk geometry in the presence of shear 
viscosity. The energy-momentum tensor for (3 = 1 is written by using (|2.17p as 



/ po _L 2r io 
I r 4/s t- ~pr 



rpfi 



V 









po 

3r 4 /3 










2m?o 



Po 
37*73 










2br, 









37473 



2CT)Q 



(3.11) 



We find that the metric components, g TT , g yy /T 2 , g xx have the following structure by solving the 
Einstein's equation recursively: 



fW(v) +r]0 h^(v)/T^ + p\v)/r^ + ..., 



(3.12) 



Note that the viscosity dependent terms exist at the order of r~ 2//3 and these are more important 
than the higher-order terms neglected in [9]. We are considering the late time region r»l, But to 
see the effects of the viscosity, we need to keep the terms at least to the order of r™ 2 / 3 . In this paper 
we consider the viscosity effects to the minimal order. 

Now we solve the Einstein's equation recursively. The power series that appear in the solution 
can be re-summed to give a compact form of the metric. After some hard work, we can get the late 
time 5d bulk geometry given b 



1 ^ 3 



-dr 1 + 1 + 



E 

t=i 



1 + £21 
1 __3_ 

1 _ eL 

1 3 



(1— 3di)7 



2ai 



(dx 



i\2 



dz 2 

> + — , 3.13 

z z 



2 One should keep in mind that we are looking for the late time geometry; the metric (|3.13[) is correct only to the 
order of 7 and the 0(7 2 ) contributions are not unambiguously determined. The representation of (|3.13|l is chosen since 
it makes the volume of the horizon finite. 
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where 



% , Po 2t7o 

and p = T - (3.14) 



' ~ p r 2 / 3 H r 4 / 3 r 2 ' 

Notice that the energy-momentum tensor ()3.11|) can NOT be written in terms of the whole p{r). It 
is an amusing surprise to see that the final metric nevertheless can be written in terms of p{r) (apart 
from the powers) in a compact form. This implies that the position of the horizon can be determined 
solely by the energy density. 

3.1 Macroscopic Parameters from Gravity Dual 

The Hawking temperature in the adiabatic approximation is given by 

T{t) = V2/(itzo(t)), (3.15) 

where 

Zo(r) = [3/pCt)] 1 /* (3.16) 
is the time-dependent position of the horizon. Using these, we obtain 

p= \ 2 N 2 c T\r), (3.17) 

by restoring the normalization of p: p — > AnGs/r^ ■ p. The entropy per unit co-moving volume is 
given by 



S = lim 



3 / 4x /1 , P^ 4 \ ( 1 " 3a ')T 



n i+ 



pz 4 
3 



T 2a,i 



1 2y/2rrl 



Notice that the entropy is completely independent of geometric parameters (i.e., a, ). Without con- 
formal invariance this result is not guaranteed. As we emphasized in p3], the value of S at r = oo, 
which cannot be determined by hydrodynamics alone, is precisely determined to be 



in terms of the initial condition po and the microscopic gauge theory parameter N c . Similarly, 

r W = ^(l-2 7W )'/*, with 7b (3-20) 
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These parameters So^Tq are precisely the quantities used in macroscopic theory (hydrodynamics) 
which should be provided by a microscopic theory like QCD. What we are showing here is that by 
considering the AdS/CFT dual of hydrodynamics, one can determine such quantities. 

Let us check consistency of (|3,18p and (12.24H . Apparent time dependence agrees in the leading 
order. In fact one can do more. The normalized entropy-creation rate from the gravity dual result 
(|3.18h is given by 

1 dS r)o ^_-7/3 



S dr poT 5 / 3 
and that from hydrodynamics result (|2.24p is 



+ 0(t~'^), (3.21) 



1 dS _ 4 t]q ^,.-7/3 



Sdr 3 ToS^t 5 / 3 
Comparing (|3.2ip and (|3.22|) . we obtain 

s - ^El = ^EL 

00 3 T 3 T 



+ 0{t-" 6 ). (3.22) 



(3.23) 



This is the consistency condition that is required to be checked|f| With use of Soo and To given in 
eqn's (|3.19p and (|3.20p we can check that the consistency condition (|3.23[) is indeed satisfied. 



4 Flow of RHIC fireball and Kasner spacetime 

So far, we have established the gravity dual of the Yang-Mills system in Kasner spacetime. Now we 
would like to suggest a relevance of our model to description of the elliptic flow in RHIC experiments. 
We have seen in section [2] that the hydrodynamic description of the three-dimensional expansion in 
the Kasner spacetime is as simple as that of the Bjorken expansion in the flat spacetime. Therefore, it 
is great if we can apply such a simple formalism to anisotropically expanding RHIC fireball. However, 
from the realistic point of view, we sacrificed the flatness of 4d spacetime for the simplicity of the fluid 
dynamics of 3 dimensional expansion. One immediate question is when and under which condition we 
can justify it. The hydrodynamics on the Kasner spacetime provides a well-approximated description 
of a three-dimensional expansion in the flat spacetime if curvature which is small enough. Notice 
also that the spacetime symmetry crucial to our problem, which is uniformity of the spacetime, is 
maintained in Kasner spacetime. 

In figure 1, we show the allowed region of the anisotropic parameters, (b, c) = (0,0) corresponds 
to the Bjorken expansion. The non-zero components of the Riemann tensor of the Kasner spacetime 
are 

R 0l0 i = (1 - ai)a t T 2a >- 2 , Rijij = a iaj T 2a ^- 2 , (4.1) 

3 In fact this is the relationship among the entropy, the energy (per unit co-moving volume) and the temperature 
obtained by thermodynamics at r = oo where the system reaches thermal equilibrium. 
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-0.5 



Figure 1: The available region of (b,c) is on the ellipse between (0,0) and (1,0). (6, c) = (0,0) is Bjorken 
point 

and the non-flatness is directly related to the distance from the Bjorken point (a, b, c) = (1,0,0) 
on the parameter space. In order to keep the deviation from the flat spacetime small, we restrict 
ourselves within the vicinity of the Bjorken point, 

a~l, 6~0, c~0, (4.2) 

that corresponds to almost central collisions in RHIC. 

In appendix, we show what approximation is necessary to reach kasner space starting from a flat 
spacetime. The conditions are: 1) the fluid is produced by almost central collision (small b), 2) we 
consider only the central part of the fluid (small X -1- ), and 3) we consider only the late time regime. 
For more detail, see appendix \M 

With these limitations in mind, let us consider how the elliptic flow can be described within our 
framework. We can choose b > c without any loss of generality. By considering the intersection of 
the plane and the unit sphere in the a,b,c space of (|2.3[) . we can see that 



a > b > -c> 0. (4.3) 

This means that one of the transverse directions must contract and the others expand so that the 
expansion is elliptical whose eccentricity grows and eventually saturate to 1. @ 

The transverse expansion has a natural interpretation as the elliptic flow, one of the most concrete 
evidence for the strong nature of the interaction. Let e be the eccentricity defined b}@ 



4 There is no contraction in the realistic RHIC QGP. The contraction in the present model is due to the conformal 
invariance which is unavoidable for TV = 4 SYM theory. Notice that this volume-preserving nature indicates our TV = 4 
SYM plasma is more "liquid-like" than the RHIC QGP. 

The e defined here contains an extra minus sign comparing to other literature such as [151 116] . We define its 
positivity by the direction of the V2 evolution. 
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where X % is the cartesian coordinate and (• • -)x = / • • • p dX2<iX^. On the other hand, V2, the 
quantity experimentally characterizing the elliptic flow, is defined by (see for example, |X5 |, 116 ^ [20]) 

4tTT = T~ ("o + 2v 2 cos(20) + 2u 4 cos(40) + •••), (4.5) 
Jy dip 2tt 

where N is the number of the partons and (p is the angular coordinate on the transverse momentum 
plane. It can be calculated from the following identification, 

2_ p2 > 



fd 2 P. ( P (- P * \ dN 



*> - s ^£_ ■ m 

where P % is the momentum of the fluid in the X 1 coordinate and d 2 Px = dP 2 dP s . To consider v-i 
in the present model, we introduce a coarse-grained (i.e., averaged over a small volume) momentum 
flow, 

P^KHx^t). (4.7) 

At each fixed time r, this can be considered as a mapping from the P l -space to the a^-spacejj Now 
i>2 can be expressed as an integral over the co-moving coordinates: 

2_ p2 s 



Jdx 2 dx 3 ( pl + % ) Pn(x) 



V2 = n ,% , (4.8) 

J dx^dx pn{x) 

where pn{x) = d ^ dx i is the particle density in the transverse space. Notice that the Jacobians cancel 
out. 

To proceed, we need an explicit expression for K(x). We first relate the coordinates of the Kasner 
spacetime and the usual Minkowski spacetime. We can identify the flat space variable X 1 for i = 2, 3 
by 

X i ~ r at x i near x i = 0, (4.9) 
which is nothing but (IA.2j) under (IA.5j) . Then the fluid momentum in the flat space is given by 

P i = p^— = paiT^xK (4.10) 
ar 

Here p = mp^ is a mass density with some proper mass parameter m. Since we treat it as a constant 
from now on, it is irrelevant in the calculation below. 

A few technical remarks are in order: 

1. One should notice that the fluid momentum P l (x) is different from the individual particle 
momentum. By replacing the momentum by fluid momentum, we expect a small deviation from 
the original v-i- So one may want to call the final expression by V2- However, this is precisely 
what we should have when we describe the system by hydrodynamics where everything is to 
be defined in the coarse-grained level. 



In the co- moving coordinate, where the fluid is at rest, we do not have any flow. So P l (x) is the coarse-grained 
momentum field in Minkowski space written as a function of co-moving coordinate. 
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1 2 3 4 5 



Figure 2: Time evolution of the elliptic flow parameters. Dotted line is the eccentricity and the solid line is 
the V2 for b — 0.1. 



2. Since our approximation is valid in the small x 1 region, the integrals in (|4,8p are now defined 
within the small x % region by introducing a cut off radius. 

The resulting V2 can be calculated in our model to yield 
for b > 0, c < 0. Comparing this with the result for the eccentricity in the small x % region 

Jib _ 2c 

' ' (4-12) 



T 2b i T 2c 



we plot the time evolution of V2 in figure 2. 



Notice that the hydrodynamics describes relatively late time regime and we do not consider the 
time region where V2 is negative. At the time of zero V2, £ is negative and its absolute value decreases 
to zero. That is, the fireball core becomes more round in the transverse space. After that initial 
period, e becomes positive and grows in the same direction of growing v%. This is qualitatively 
the same behaviour as those given in |15t [T6] and [20] . Notice that there is only one independent 
parameter b that can be used to parameterize the initial eccentricity. 



5 Discussion 



In this paper, we extended the Bjorken hydrodynamics to the case of anisotropic three-dimensional 
expansion. The Ricci flat condition suggested by the conformal invariance imposes the condition 
on the anisotropy parameters. As a consequence, the four-dimensional boundary metric becomes 
precisely that of Kasner universe. Our hydrodynamic equation of motion is independent of the 
anisotropic parameters, so that it is the same as that obtained by Bjorken. 

Although the Kasner spacetime is a curved spacetime, we found that it gives a well-approximated 
local rest frame of an anisotropically expanding fluid on the flat spacetime. We obtained the eccen- 
tricity and v 2 of the elliptic flow of the fluid. The expansion in this setting has a deviation from that 
of realistic RHIC fireballs, since one of the transverse direction contracts as a consequence of the 
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Kasner condition. This deviation is essentially due to the conformal invariance. Whether one can 
lessen the condition of the conformal invariance is an important issue, which we will treat in other 
publication. 

We also extended the "falling horizon solution" obtained in (9j[H] to the case of three-dimensionally 
expanding fluid. It is very important to figure out how to use the bulk metric apart from reading off 
the horizon location, which gives the cooling rate and the entropy creation rate. It is also interest- 
ing to workout the details of the Hawking evaporation and Wilson line calculations for the external 
quark-antiquark in our time-dependent metric. We will come back to these issues in later publication. 
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A Relationship between Kasner spacetime and our local rest frame 

We clarify the relationship between our local rest frame and the Kasner spacetime in the following 
strategy. We start from the cartesian coordinate and make a coordinate transformation to define our 
local rest frame. The transformation is essentially what we have employed at (14. 9h , At this stage, we 
are still on the flat spacetime but the local rest frame is that of the anisotropically expanding fluid 
with elliptic flow. Next, we see that the local rest frame is equivalent to the Kasner spacetime under 
a well-controled approximation which is specified at the end of this section. 

Let us start from the cartesian coordinate, 

ds 2 = -{dX°f + (dX 1 ) 2 + (dX 2 ) 2 + (dX 3 ) 2 , (A.l) 

and perform the following coordinate transformation 

(X^X^X^X 3 ) = (Tcoshx 1 ,rsinhx 1 ,r fe x 2 ,r c x 3 ) (A.2) 

to define our local rest frame (r, X j X <j X 3 ). Notice that the (X^X 1 ) part has been "boosted" with 
rapidity x 1 and the (X 2 ,X 3 ) part has been transformed by (14. 9|) . The resulting metric on the local 
rest frame is given by 

ds 2 = -dT 2 + T 2a (T 2 - 2a )(dx 1 ) 2 +T 2b ^dx 2 + b^dT^ + T 2C (dX 3 + CydT^J . (A.3) 

The parameters a, b, c are free so far. 
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Next, let us impose the Kasner condition (|2.3p among a, b, c and assume b <C 1. a and c is then 
given by a = 1 — 36 2 and c = — 6 + 36 2 approximately. Now (|A.3|) is written as 



ds 2 = -dr 2 + r 2a (l + 6b 2 log r + • • •) (dx 1 )' 

+T 2b (dx 2 + b—dr] + r 2c I dr 4 + ( -/> 3/r + • • • ) — ,■/- I . (A.I) 



The above metric agrees with the Kasner metric (|2.2p under the following approximation: 

I X ^ I 1 3> ^ I 

J — l ~ J — I ~ b < 1, with 61ogr ~ 0(1). (A.5) 
r r 

We ignore 0(b 2 ) contributions while we keep felogr. This is a mixture of small b approximation (we 
consider almost central collision), small (x 2 ,x 3 ) approximation (we consider only within the central 
part of the QGP) and the late time approximation. The condition Mogr ~ 0(1) is necessary for us 
to get the nontrivial results. 
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